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ABSTRACT 

An approximate statistical description of the formation and evolution of structure of 
the universe based on the Zel'dovich theory of gravitational instability is proposed. It 
is found that the evolution of DM structure shows features of self-similarity and the 
main structure characteristics can be expressed through the parameters of initial power 
spectrum and cosmological model. For the CDM-like power spectrum and suitable 
parameters of the cosmological model the effective matter compression reaches the 
observed scales Rwaii ^'^0 - 25/i~"'^Mpc with the typical mean separation of wall-like 
elements Dslss ~50 - 70/i~^Mpc. This description can be directly applied to the 
deep pencil beam galactic surveys and absorption spectra of quasars. For larger 3D 
catalogs and simulations it can be applied to results obtained with the core-sampling 
analysis. 

It is shown that the interaction of large and small scale perturbations modulates 
the creation rate of early Zel'dovich pancakes and generates bias on the SLSS scale. 
For suitable parameters of the cosmological model and reheating process this bias can 
essentially improve the characteristics of simulated structure of the universe. 

The models with 0.3 < flm < 0.5 give the best description of the observed struc- 
ture parameters. The influence of low mass "warm" dark matter particles, such as a 
massive neutrino, will extend the acceptable range of flm and h. 

Key words: cosmology: large-scale structure of the Universe — galaxies: clusters: 
general - theory. 



1 INTRODUCTION 

Over the past decade the large maps of the spatial galaxy 
distribution have been prepared and the unexpectedly com- 
plicated character of this distribution was established. The 
structure predicted by the Zel'dovich theory of gravitational 
instability (Zel'dovich 1970, 1978) was found already in the 
first wedge diagrams (Gregory & Thompson 1978) and 
now the Large Scale Structure (LSS) is seen in many obser- 
vational catalogs, such as the CfA (de Lapparent, Geller & 
Huchra 1987; Ramella, Geller & Huchra 1992), the SRSS (da 
Costa et al. 1988) and in the Las Campanas Redshift Sur- 
vey (Shectman et al. 1996, hereafter LCRS). The observed 
high concentration of galaxies within the wall-like structure 
elements such as the Great Attractor (Dressier et al. 1987), 
and the Great Wall (de Lapparent, Geller & Huchra 1988) 
and the existence of extended under dense regions similar 
to the Great Void (Kirshner et al. 1983) put in the forefront 
the investigation of the Super Large Scale Structure (SLSS). 



Now the SLSS is also found in many deep pencil beam red- 
shift surveys (Broadhurst et al. 1990; Willmer et al. 1994; 
Buryak et al. 1994; Bellanger & de Lapparent 1995; Cohen 
et al. 1996) as a rich galaxy clumps with the typical sepa- 
rations in the range of (60 - 120)/i"^Mpc. Here h=//o/100 
km/s/Mpc is the dimensionless Hubble constant. 

Further progress in the statistical description of the LSS 
& SLSS has been reached with the core-sampling method 
(Buryak et al. 1994) and the Minimal Spanning Tree tech- 
nique (Barrow, Bhavsar & Sonoda 1985). Recent analy- 
sis of the LCRS performed by Doroshkevich et al. (1996 
& 1997b, hereafter LCRSl & LCRS2, Doroshkevich et al. 
1998) revealed some statistical parameters of the wall-like 
SLSS component such as their typical separation, Dslss ~ 
50 — 60/i^^Mpc, and the fraction of galaxies accumulated by 
the SLSS, which can reach ~ 50%. The same analysis indi- 
cates that formation of richer walls can be roughly described 
as an asymmetric 2D collapse of regions with a typical size 
Rwaii ~ 20 - 25/i~^Mpc that is about half of their typical 
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separation. The analysis of Durham/UKST redshift survey 
confirms tliese results (Doroshkevich et al., f999). Earlier 
similax scales, in the range of (50 -100) ft" ^Mpc. were found 
only for spatial distribution of clusters of galaxies (sec, e.g., 
Bahcall 1988; Einasto ct al. 1994) and for a few superclusters 
of galaxies (see, e.g., Oort 1983 a,b). 

Evolution of structure was discussed and simulated 
many times (see, e.g., Sahni et al. 1994; Doroshkevich et 
al. 1997a, hereafter DFGMM; for references, Sahni & Coles 
1995). However, SLSS in the dark matter (DM) distribution 
similar to that seen in the LCRS was found only recently 
in a few simulations with the CDM-like power spectrum 
and nmh= 0.2 - 0.3, (Cole, Weinberg, Frcnk & Ratra 1997; 
Doroshkevich, Miiller, Retzalf & Turchaninov 1999, here- 
after DMRT). Hence, for suitable cosmological models the 
evolution of small initial perturbations results in the SLSS 
formation. 

In this paper we present an approximate statistical de- 
scription of the process of DM structure formation based on 
the nonlinear Zel'dovich theory. The potential of this ap- 
proach is limited as the successive consideration of mutual 
interactions of the small and large scale perturbations be- 
comes more and more cumbersome. In spite of this it allows 
us to obtain some interesting results. Thus, it is shown that 
formation of both LSS and SLSS is a joint process possess- 
ing some features of self - similarity. The main observed 
characteristics of LSS and SLSS are expressed through the 
structure functions of power spectrum and through the typ- 
ical scales, set by the power spectrum, the time scale, set by 
the amplitude of perturbation, and the main parameters of 
cosmological model. One of the most interesting such charac- 
teristics is the dynamical scale of the nonlinearity defined as 
the scale of essential DM concentration within high density 
walls. We show that for the CDM transfer function (Bardeen 
et al. 1986, hereafter BBKS) and Harrison - Zel'dovich pri- 
mordial power spectrum and for cosmological models with 
lower matter density this scale of nonlinearity reaches 20 
- 30/i~^Mpc that is comparable with typical scales of the 
observed SLSS elements. 

Simulations (DMRT) show that even in cosmological 
models with a low matter density the sinmlated velocity 
dispersion within the SLSS elements reaches 400 - 700km/s 
along each principal aocis; this exceeds the observed value by 
a factor of ~ 1.5 - 2. Such a large and isotropic velocity 
dispersion is caused by the disruption of the walls into high 
density clouds. For smaller matter density of the universe 
this dispersion decreases but together with the fraction of 
matter accumulated by the walls. This means that other 
factors as, for example, the large scale bias in the spatial 
galaxy distribution relative to the more homogeneous dis- 
tribution of DM and baryoiis could be essential for the suc- 
cessful reproduction of the observed SLSS. Such large scale 
bias caused by the interaction of small and large scale per- 
turbations was discussed by Dekcl & Silk (1986), and Dekel 
& Rees (1987), and estimated by Demianski & Doroshke- 
vich (1999, hereafter Paper I). 

The interaction of small and large scale perturbations 
is important during all evolutionary stages. Thus, even dur- 
ing early evolutionary periods the large scale perturbations 
modulate the rate of pancake formation. This modulation 
is seen as an acceleration of the pancake formation within 
deeper potential wells which later are transformed into the 



wall-like SLSS elements (Buryak et al. 1992; Demianski & 
Doroshkevich 1997; Paper I). Suppression of pancake forma- 
tion near the peaks of gravitational potential noted by Sahni 
et al. (1994) is another manifestation of such interactions. 
During all evolutionary stages these interactions result in 
the successive merging of individual pancakes. The acceler- 
ation of pancake disruption, caused by compression of mat- 
ter within walls, can also be attributed to this interaction. 
Now it is observed as a high velocity dispersion in simulated 
SLSS and as differences between the expected and measured 
mass functions. It was found to be essential even for pan- 
cakes formed at high redshifts (Miralda-Escude et al. 1996). 
The possible correlation of galaxy morphology with large 
scale perturbations was discussed by Evrard et al. (1990). 
All these manifestations of small and large scale interaction 
are important for the correct comparison and interpretation 
of simulated and observed matter distribution. 

Now the modulation of spatial distribution of pancakes 
formed at high redshifts z > 4 can be seen as the large scale 
bias in the galaxy and DM spatial distribution. This bias 
can be generated by the combined action of large scale per- 
turbations and reheating of baryonic component of the uni- 
verse (sec, e.g., Dekcl & Silk 1986; Dekcl & Rccs 1987). The 
reheating was discussed many times during the last thirty 
years in various aspects (see, e.g., Sunyaev & Zel'dovich 
1972; White & Rees 1978; Shapiro, Giroux & Babul 1994). 
Effects of reheating on the process of galaxy formation were 
discussed as well (see, e.g.. Babul & White 1991; Efstathiou 
1992; Quinn, Katz & Efstathiou 1996). It is also known 
that under reasonable assumptions about the possible en- 
ergy sources reheating can occur for relatively small range 
of redshifts z ~ 5 — 10 (sec, e.g., Tcgmark et al. 1997; Baltz 
et al. 1998). If essential concentration of baryons in high 
density clouds is reached at the same redshifts, the reheat- 
ing can help to generate bias (Demianski & Doroshkevich 
1997; Paper I). In this case further formation of high den- 
sity baryonic clouds will be significantly depressed, due to 
reheating, within extended regions observed today as under 
dense regions between richer walls. Our estimates show that 
this spatial modulation of the luminous matter distribution 
may be essential for the interpretation of observations. 

Numerical simulations are now the best way to repro- 
duce and to study the joint action of all the pertinent fac- 
tors together and to obtain more representative descrip- 
tion of the process of structure formation. Essential progress 
achieved recently both in the simulations and study of DM 
and 'galaxy' distributions (Governato et al. 1998; Jenkins 
et al. 1998; Doroshkevich, Fong & Makarova 1998; DMRT; 
Cole, Hatton, Weinberg & Frenk 1998) allows us to follow 
the structure evolution in a wide range of redshifts and to re- 
veal differences between DM and galaxy distribution. Com- 
parison of these results with observations and an approxi- 
mate theoretical description stimulates further progress in 
our understanding of evolution of the universe. 

This paper is orgaiuzed as follows: In Sect. 2 main nota- 
tions arc introduced. In Sec. 3 and 4 the distribution func- 
tions of DM pancakes are derived and the interaction of 
small and large scale perturbations is described that allows 
us to obtain in Sees. 5 and 6 the statistical characteristics of 
DM structure. In Sec. 7 the large scale bias is discussed and 
in Sec. 8 the dynamical characteristics of walls are found. In 
Sec. 9 the theoretical estimates are compared with the avail- 
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able observational and simulated data. We conclude with 
Sec. 10 where a short discussion of the main results is pre- 
sented. Some technical details are given in Appendixes I - 
IV. 



2 STATISTICAL PARAMETERS OF 

PERTURBATIONS: VARIANCES AND 
TYPICAL SCALES 

The simplest characteristics of perturbations are the vari- 
ances of density and velocity perturbations. For a more do- 
tailed statistical description of the structure evolution it is 
necessary to use also the structure functions. They were in- 
troduced in Paper I and are briefly described in this Section 
and Appendix I. Here we consider only the SCDM-like power 
spectrum but the same approach can be applied for other 
spectra as well. 

Our analysis is based on the Zel'dovich theory which 
links the Eulorian, r^, and the Lagrangian, qi, coordinates 
of fluid elements (particles) by the expression 

n = {l + z)-%-Biz)S,ici)], (2.1) 

where z denotes the redshift, B{z) describes growth of per- 
turbations in the linear theory, and the potential vector 
Si{<l) = d(j)/dqi characterizes the spatial distribution of per- 
turbations. The Lagrangian coordinates of a particle, qi, are 
its unperturbed comoving coordinates. 

For the flat universe with Q,m = 1, f2m > 0.1, the 
function B{z) can be approximated with a precision better 
then 10% by the expression (Paper I) 

l-nm+2.2Q.m(l + zf 



^ ■ 1 + 1.2^^ 

and for an open universe with 0.1 < n„ 



(2.2) 



< 1, Q,A = as 

(2.3) 
= both 



1 -I- 1.5^™ 

(Zel'dovich & Novikov 1983). For 0,^ = 
expressions give B~^{z) = \ + z. 

The main characteristics of the perturbations are the 
variances of density (7^, displacement crj, and components 
of the deformation tensor g% 



= ^ y p{k)dk, 



: 5ct|, 



^ j P{k)k dk, 
(2.4) 

where p{k) is the power spectrum, and k is the comoving 
wave number. 

The power spectrum determines also two amplitude in- 
dependent typical scales, which allow us to describe the pro- 
cess of structure formation and can be, possibly, estimated 
from the observed galaxy distribution. For the Harrison - 
Zel'dovich primordial power spectrum these scales, lo and 
Ic, are defined as 



Jo 



kr{k/ko)dk, 
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where T^{x) is the transfer function and ko = Qrnh?yipc ~^ . 
For the CDM transfer function (BBKS) the scale lo and the 
typical masses of DM and baryonic components associated 
with the scale lo are 



lo « 6.6(nmft)' 
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For the SCDM power spectrum the Up is divergent (log- 
arithmically) at large k and the standard CDM transfer 
function should be truncated by introduction of an appro- 
priate mass of DM particle, Mdm- This restriction allows, 
however, a significant interval for a possible mass of DM 
particles - up to lO^** - 10^" eV. In this paper wc will use 
the ratio Ic/lo as a measure of the mass of DM particles. For 
0.5 keV <Mdm< lO" keV we have 0.2 > U/lo > 0.01 and 
for Mdm > 5keV the following approximate expression can 
be used 

Wlo^^^, (2.7) 



M,= (ff Mo 



Vln Mdm 

^ 0.6 ■ 10'°AfQ 
(lnMDM)-V^0?„/t4' 



where Mdm is the mass of DM particle in keV. It turns 
out that the final results only weakly depend on the mass of 
DM particles. For the Harrison-Zol'dovich primordial power 
spectrum the dependence of the moments and typical scales 
on the mass of DM particles was discussed in Paper I. More 
cumbersome correlation and structure functions of pertur- 
bations were also discussed in Paper I and, partly, are pre- 
sented in the Appendix I. 

The amplitude of power spectrum can be taken from 
the measured anisotropy of the relic radiation (Stompor et 
al. 1995; Bunn & White 1997; Gorski et al. 1998). It can be 
fitted by 



; 31.2^ 



a{flm, Ha = 1 — fir, 

a{flrn, = 0) — Q 



,0.215-0.06infJ„ 



0.65-0. 19inn„ 



where Tq is the amplitude of quadrupole component of 
anisotropy of the relic radiation. The time scale of the struc- 
ture evolution is defined by the function 



r{z) = ToB{z), 



(2.9) 
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, ^A = 0, < 1. 



Further on, as a rule, the dimensionless variables will be 
used. We will use lo as the unit of length, and < p > lo as 
the unit of surface density. This means also that such dimen- 
sionless characteristics of a pancake as the size of collapsed 
slab and resulting surface density of a pancake are identical. 
Below we will use both terms as well as the term 'mass', m, 
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to characterize the surface density reached during formation 
of a pancake. The gravitational potential and displacement 
are measured in units of aslo/V^ and as/VS, respectively. 



3 STATISTICAL CHARACTERISTICS OF 
PANCAKES 

In this section the mass function of Zel'dovich pancaJtes and 
its time evolution is given. This can be done using the main 
equation of Zel'dovich theory (2.1). The mass of compressed 

matter is measured by the Lagrangian size of compressed 
slab, q, or by the surface density of pancake, < p > g. As 
it was noted above, both measures are identical in dimen- 
sionloss notation. So, wc will use both terms 'size', q, and 
'mass', rn, of a pancake to characterize its surface density. 

Here we do not consider the transversal characteristics 
of structure elements and cannot discriminate, for example, 
the central part of a poorer pancake and periphery of richer 
pancake, if they have the same surface density or mass, m. 
In this sense our approach gives characteristics similar to 
that obtained with the core-sampling analysis, pencil beam 
observations or the distribution of absorption lines in spectra 
of QSOs. 

3.1 The pancake formation 

According to the relation (2.1) when two particles with dif- 
ferent Lagrangian coordinates qi and q2 meet at the same 
Eulerian point r a pancake with the surface mass density 
< p > |qi— q2| forms. Here we assume that aZZ particles situ- 
ated between these two boundary particles are also incorpo- 
rated into the same pancake. This assumption is also made 
in the adhesion approach (see, e.g., Shandarin & Zel'dovich 
1989). Formally, this condition can be written as 

qi2 = qi - q2 = r • [S(qi) - S(q2)]. (3.1) 

This means that the pancake formation process can be char- 
acterized by the scalar random function 



Q(gi2) = ^ • [S(qi) - S(q2) 
912 



(3.2a) 



under the condition 



qi2 X [S(qi) - S(q2)] = 0. 

In a coordinate system with the first axis oriented along 
the vector qi2 the condition (3.1) can be rewritten more 
explicitly as 

Q(gi2) = A5i = gi2/T, A52(gi2) = AS'3(gi2) = 0, (3.26) 

AS'i = Si(qi)-&(q2), 

and formation of a pancake is described by the ID matter 
flow between points qi — qi ■ qi2/'7i2 and 172 = q2 ■ qi2/gi2. 

As is shown in Appendix II, under the assumption of 
Gaussian distribution of perturbations and neglecting the 
(weak) correlations between matter motion in orthogonal di- 
rections, the probability of a pancake formation, for a given 
q prior to the 'time' r or for a given r with a size larger 
then q, is identical to the probability to have ASi > q/r 
(Q(9i2)/gi2 > 1/t), so 

(MY' 



Wcr{> q,T) 



1 + erf 



\V2r J 



= , ^ 

V2[l-Gi2(g)] 

where (see Appendix I &i II) 

M(g) « , q<qo, 

H{q) ~ V9/2, go < 9 < 1, m(«) « g/V2, Q > 1, (3.4) 
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As is seen from (3.1), and (3.3) the pancake characteristics 
at the moment r are expressed through the function fjb{q) 
related to the initial power spectrum. This is an essential 

&>aturc of the Zcl'dovich theory allowing one to obtain ap- 
proximate analytical description of pancake properties. 

Equation (3.3) demonstrates that ~ | of matter with 
Q{<li2) > is compressed at least in one direction and for 
~ i of matter 3D expansion takes place what reflects the 
symmetry of initial distribution of the displacements, A/Sj. 
At first glance it seems that this conclusion is in contrar 
diction with known results of Zel'dovich theory. Actually, 
the fraction of mass, /dm, accumulated by pancakes can 
be estimated using the PDF for the largest principal value 
of deformation tensor, dW{\\), (Zel'dovich 1970; Doroshke- 
vich 1970; Doroshkevich & Shandarin 1979; Shandarin & 
Zel'dovich 1989; Paper I). This PDF can be approximated 
(with a precision ~10%) by the Gaussian function with 
< Ai >= ?,gd/\/2^, a'i = (13/6 - 4.5/7r)a|). As in the 
Zel'dovich theory pancake formation is described by the re- 
lation B{z)Xi = 1 (which follows directly from (2.1) ), so for 
the fraction of compressed matter with Ai > 1/B we have 



foM w -^erfc 



V2t 



(3.5) 



(3.3) 



and for r 3> Ic/lo, foM ^1 (in the Zel'dovich theory 
0.92< foM <1). This shows that already during the early 
period of nonlinear evolution, at r w Ic/lo ^ 1, large frac- 
tion of matter foM > 0.9 is compressed into low mass pan- 
cakes. But for the CDM-like power spectrum the description 
of matter compression through the deformation tensor is ap- 
propriate only at small scales q < qo whereas for q ^ qo, 
the correlations between matter flow in orthogonal direc- 
tions rapidly decrease what can be seen directly from the 
expressions for the structure functions given in Appendix I. 

At larger scales we have to use the more cumbersome 
description discussed above and the estimate (3.3) shows 
that at such scales the efficiency of matter integration into 
structure elements is only ~ 0.875 (more accurate estimates 
taking into account the correlation of displacements lowers 
this value to 0.79). This limit is reached already at small r, 
for g -C 1, that means strong matter concentration within 
small structure elements. Further evolution does not change 
this limit and only redistributes - duo to sequential merging 
- the compressed matter to more and more massive struc- 
ture elements. Thus, the approximate estimates show that ~ 
21% of matter is subjected to 3D compression, ~ 21% to 3D 
expansion, ~ 29% of matter is subjected to 2D compression 
and can be accumulated by filaments and ~ 29% is subjected 
to ID compression and remains in pancakes. The difference 
between estimates (3.3) and (3.5) shows that ~ 15 - 20% 
of matter incorporated in small clouds, with M < Mq, is 
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not accumulated by larger pancakes, with M > Mq, and re- 
mains distributed between those pancakes. These estimates 
can be changed because the Zel'dovich approximation be- 
comes invalid when strong matter compression is reached 
during pancakes formation. 



When Ic <C lo it is described by the following (approximate) 
expression: 



3.2 The characteristics of pancake formation 

The probability distribution function (PDF) for pancakes 
formed at the moment r can be found from (3.3) as 

7 dq 7V27rr dq \ V2t / 

(3.6) 

$(a;) = e-^'[l-|-er/(a;)]^ 

Simple analysis based on relations (3.3), (3.6) shows that 
the maximum of the PDF (3.6) is reached at g « 1.5go <d. 

The mean mass (surface density) of formed pancakes is 
described by the expression 



< m(r) >= / qNcrdq = | / 
Jo Jo 



qNcrdq =^ j Wcrdq, (3.7) 



lcn{> q) 



1.35 m(9)^ ( 



q 



(3.10) 



where < Hr >= 0.5 describes the random orientation of a 
pancake and the chosen straight line (see also Appendix II). 

To characterize the 3D distribution of higher peaks of 
the function Q{q) and the percolation process, in the course 
of which separate pancakes are integrated into the joint 
structure, the standard technique using the Euler charac- 
teristics can be applied. It is described in details by BBKS; 
Gott et al. (1989); Coles et al. (1996); Seto et al. (1997). 

In our case the anisotropy of the pancake compres- 
sion manifests itself as an anisotropy of correlation func- 
tions of displacements and their derivatives what re- 
sults in the appearance of different Euler characteristics, 
"31 (g, r ) , n32 (g, r) , nss (?, r) : 



nai 



3 < Q22Q33 - Qlz > 
327r2 



< IQil > * 



V2t 



< m(T) >Ri 4t , r < 1, < m(T) r, r > 1, 
and the mass distribution is characterized by the function 



l^J^T <m> dq \ 
The rate of formation of pancakes with mass q is 



. ^ 8 dWcr , . 



" \V2r 



(3.8) 



(3.9) 



and the maximal rate is reached for /x(g) « t. 

These relations take into account the growth of pan- 
cakes size due to accretion of matter and form a reliable 
basis for further theoretical considerations of various char- 
acteristics of structure elements. However they cannot be 
applied directly to the observed and simulated matter dis- 
tribution as they do not consider the possibility of merging 
of pancakes. In the process of merging, earlier formed pan- 
cakes are accumulated by larger ones. A more refined tech- 
nique taking into account the pancake interaction should be 
used in order to describe the evolution of structure elements 
that are the pancakes that survived the merging process. 
The same problem appears in the Press-Schechter theory 
and it can be solved by introduction of the survival prob- 
ability (Peacock & Heavens 1990; Bond et al. 1991). This 
problem will be discussed below. 



3.3 Spatial chciracteristics of pancake distribution 

Using the standard technique (see, e.g., BBKS) we can also 

find the mean cumulative comoving linear number density 
of pancakes, n(> q), which characterizes the pancake distri- 
bution along a random straight line 



n{> q)dr = 



3 l^jq) 
4V27r q 



< VQ • dr > $ ( ^ 



Qi 



dQ 

dqi ' 



Qi. 



dqidqj 



(3.11) 



where variances an, adi, 0-^2, fds are defined by (II. 4) and 
(II.6) and expressions for 7132 = 7133 can be found by cyclic 
permutation of indices of derivatives Qi&iQij. 

Functions 713-2 & '133 characterize the percolation pro- 
cess in the direction qi2, and nsi characterizes this pro- 
cess in transversal directions. For the most interesting case 
go <C 1, go < g we have 



l3 

«0'^32 



^0^33 



;3 

10^31 



M 1 + 2g 
"°l7^-3(TT^ 



^ no 



(3.12) 



no 



3 
47r2 



2 I ggo 



3/2 



\V2t 



The function nsi is similar to the standard expression for 
an isotropic Gaussian field (BBKS). 

The first zero of Euler characteristic describes ap- 
proximately the percolation when separate higher peaks 
are incorporated into a larger (in the limiting case - in- 
finite) structure clement (Tomita 1990; Mccke and Wag- 
ner 1991). The expressions (3.12) show that in the direc- 
tions orthogonal to qi2 the percolation takes place at the 
'moment' r = /i(g) whereas along qi2 it occurs later, at 
r — /i(g)\/3(l + g)/(l + 2g). This fact favors two step per- 
colation and percolation along the pancake surface takes 
place first. For a given g and r the ratio n{> g)/n3i and 
n(> g)/n32 measure the mean surface of formed pancakes. 

Expressions (3.7) & (3.12) show that at small r perco- 
lation takes place for pancakes with m ^ < m > w 4t^ . 
This means that only the highest peaks with m > < m > 
, > T can be considered as discrete objects. The number 
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density of such peaks can be obtained from (3.12), for /x/r > 
1, as 



3 M 

lonpk{> q,T) = no^. 



(3.13) 



Later, when objects with m ~ 1 are formed and an essential 
fraction of matter is concentrated in larger pancakes with 
m ~ < m > > T, the percolation takes place through 
relatively low mass pancakes with m < < m>. This ideal 
picture will be however strongly distorted by the matter 
compression or expansion in transversal directions and by 
pancake disruption due to the gravitational instability. 



3.4 Pancake interaction 

An important characteristic of pancake interaction is the 
two points PDF which gives the conditional probability to 
form a pancake with the size q at the moment t at the 
distance T>sep from a pancake with the size Di formed at 
the moment Ti. Here we will restrict our consideration to 
the simplest case when all three vectors, q, Di,&:Dsep are 
aligned along the same straight line that is the most interest- 
ing case for the small separations. Much more cumbersome 
case of arbitrary axrangement and/or orientation of these 
vectors will be considered separately. 

In this simplest case the conditional probability to have 
AS > q/r is 



(> q, t; ^,D,^pj = 0.5 • erfc{g2/V2), 



(3.14) 



9{x,y,rs) 



92 = g\p{q)/T,n{Di)/Ti,rs]. 



The function /i(g) is given by (3.3) and the coefficient of 
correlation, Vsiq, -Di, Dsep), describing the interaction of two 
pancakes, can be found in Appendix III (III.l). The condi- 
tional distribution function of such pancakes can be written 
as 



\T Ti ) y 7T dq 



-0.59^ 



(3.15) 



These relations characterize the parameters of pancakes 
formed at the moment t under the condition of a pancake 
formation at the moment n with the size Di and the sepa- 
ration Dsep- 

The basic relation (3.1) implies that two pancakes with 
sizes Di and D-z and a separation \Dsep\ < 0.5(Z)i + D2) 
merge together and form a single pancake. For larger separa- 
tions merging of pancakes can also be considered in the same 
manner as before, but using the Euler position of formed 
pancake 



qcent - B{z) 



A01 



qi2 



(3.16) 



where qcmt = 0.5(gi +92) is the Lagrangian position of cen- 
tral point of the pancake and A4>i2 = 4'{<li) ~ ^(92). As it is 
apparent from (3.16), rpan depends on the potential differ- 
ence, [(t>{qi) — (f){q2)]/qi2, and therefore the pancake evolution 
depends on the potential distribution. 

As discussed in Paper I around each pancake there is 
a potential well with a typical size exceeding the mass m 



of collapsed pancake by a factor of 1.5 - 2. The matter in- 
fall into these potential wells stimulates merging of pancakes 
formed with a small separation. For larger separation the in- 
fluence of this factor decreases progressively and the process 
of merging becomes (almost) random. These results are con- 
sistent with predictions of the adhesion model (Shandarin & 
Zel'dovich 1989). 

The conditional probability that two pancakes with 
sizes Di and D2 whose separation is |Dsep| > 0.5(-Di + D2) 
merge at the moment r is defined by the condition 

and instead of eq. (3.3) we have for the probability of merg- 
ing 

WmerADl, D2, Dsep) = \erfc ( xiDl, D^Dsep) ^ ^3^^^ 

The function xiDi, D2, Dsep) is given by (III. 6). We can ex- 
tend applicability of the formula (3.17) for small separations 
by requiring that 

Wmerg{Dl,D2,Dsep) = 1, \Dsep\ < 0.5(Dl -^Da). (3.18) 



4 PANCAKE EVOLUTION AND FORMATION 
OF FILAMENTS 

Similar technique can also be used to estimate the transver- 
sal size, the pancakes compression and/or expansion in 
transversal directions, and other properties of pancakes. As 
we are interested in the formation of structure elements with 
typical sizes g S> go the local description through the de- 
formation tensor cannot be used and the imposed condi- 
tions make even approximate description of pancake evolu- 
tion quite cumbersome (see, e.g., Kofman et al. 1994). The 
general tendencies and rough characteristics of this evolu- 
tion can only be outlined. Thus, for example, we can esti- 
mate the matter fraction compressed within filamentary-like 
elements and high density clumps as of about 50% whereas 
only ~ 29% of matter is subjected to ID compression. This 
estimate implies that larger pancakes could also incorporate 
an essential fraction of filaments and clumps. 

The formation of filaments as well as pancakes disrup- 
tion are stimulated by the growth of density in the course of 
pancakes compression, and therefore, probably, during early 
evolutionary stages, filaments represent the most conspicu- 
ous elements of the structure. As was discussed in Sec. 3.3, 
filaments merge to form a joint network. Later, when larger 
pancakes are formed, the evolution of pancakes becomes 
slower and disruption of pancakes dominates. These expec- 
tations are consistent with the observed and simulated mat- 
ter distribution. Thus, the conspicuous filaments are seen 
even at 2: = 3 (see, e.g., Governato et al. 1998; Jenkins et al. 
1998) whereas disrupted walls dominate at small redshifts 
(LCRSl, LCRS2, DMRT). 

4.1 The characteristics of filament formation 

Some approximate characteristics of filament distribution 
can be obtained by considering the formation of filaments as 
a sequential matter compression along two principal direc- 
tions. Such two step compression results in formation of high 
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density "ridge" surrounded by a lower density anisotropic 
halo. In a coordinate system with the first and second axes 
oriented along the directions of maximal and intermediate 
compressions this process can be approximately described 
by two equations similar to (3.2): 



<3(9i2) = qii/r, Qiyii) = yn/rf. 



(4.1) 



Here vectors qi2 and yi2 and functions Q{qi2) & Q{yi2) 
describe the deformation along the first and second coordi- 
nate axes respectively and additional conditions introduced 
in (3.2) are assumed to be fulfilled. As was discussed above 
the matter compression along the second axis is accelerated 
by the pancake formation and the function Tf{z) differs from 
that given by (2.9). 

Bearing in mind these restrictions we will approxi- 
mately characterize the probability of filament formation, 
Wlr, for given q & y prior to the 'time' r & tj or for 
given r & r/ with sizes larger then q &i y, respectively, as 
the probability to have Q{qi2)/qi2 > Q{yi2)/yi2 > I/t/, 
Q(9i2)/gi2 > l/r : 



Wi{> q,T;>y,Tf) = 1-1 



3 . / Kq) 



(4.2) 



The PDF for the filaments, N[^, can be found from (4.2) as 

jU)i _ 8 dn{q)dii{y) 



N^J>{q,T;y,Tf) 



STTTTf dq dy 



■ exp 



Mg) _ My) 

" 2r2 2r] 



1 + erf 



Ky) \ 

V2TfJ 



<q<oo, <y <oo (4.3) 



More adequate characteristic of filaments is the linear 
mass density, ruf = irqy/A, defined as a mass per unit length 
of filament. Both in observed and simulated catalogues fila- 
ments usually form a less massive part of structure elements 
with mf <Cl. For such filaments g ^ 1, ^(s) ^ \/9/2, 2/^1, 
m(2/) and the PDF N^j.{mf) can be found from (4.3) 

as 



1 



27r < m/ > VI 



dx _ v2 
— e 

X 



(x+l/x) 

(4.4) 



ruf 



< ruf > 



< ruf >= Anr Tf. 



The surface density of spatial distribution of filaments 
can be obtained using the same technique as in Sec. 3.3 . 



instead of expressions (3.2) and (3.3). The PDF for such 
pancakes is 



<mp >=4t (J < 1. 

This PDF differs from (3.6) by the form of the function $(a;) 
what decreases the PDF for larger /i/r. 

As before, expression (4.5) characterizes pancakes by 
the surface mass density of collapsed matter, q. However the 
pancake surface mass density varies with time after pancake 
formation due to transversal compression or expansion that 
results, in particular, in formation of filaments. Even if these 
transversal motions do not lead to such dramatic results 
they can change drastically the observed surface density of 
pancakes. So, the current surface mass density nip = q/sp, 
where Sp describes the variation of pancake surface caused 
by transversal motions, is a more adequate characteristic of 
pancakes. 

As was discussed above this period of pancake evolution 
is not adequately described by the Zcl'dovich theory and our 
results become unreliable. To obtain qualitative characteris- 
tics of influence of these factors we can, for example, describe 
the variation of pancake's surface as 

Sp OC (1 - r*Q(?/12)/j/12)(l - t''Q{zi2)/zi2). 

Here vectors yrz and Zi2, and functions Q(yi2) & Q(zi2), de- 
scribe the deformation along the second and third coordinate 
axes respectively. The function r* differs from that given 
by (2.9) and it can depend on transversal motions. Even 
so rough consideration shows that the exponential term in 
(4.5) is eroded, and the resulting PDF becomes power-like: 



N(mp <C 1) OC Trip N{mp S> 1) oc rUp 



(4.6) 



The pancake disruption accelerates this erosion as well and 
makes the PDF more complicated. 

This discussion shows that slowly evolving pancakes 
with slower transversal motions can be separated into a 
special subpopulation for which the surface density changes 
slowly, nip m q, and the PDF (4.5) correctly describes the 
pancake distribution during the essential period of evolution. 
This subpopulation is singled out by conditions 

\Q{zi2)/zi2\ < \Q{yi2)/yi2\ < e/r, e < 1 

and the probability of existence of such pancakes is pro- 
portional to T~^. This factor describes the disappearance of 
such pancakes in the course of evolution. This subpopulation 
can be however quite rich (see., discussion in Sec. 9). 



4.2 The pancake evolution 



Separation of filaments into a special class of structure ele- 
ments requires also a redefinition of population of 'true' pan- 
cakes. Now the probability of 'true' pancake formation de- 
fined by conditions Q(gi2)/gi2 > l/r, Q(t/i2)/2/i2 < 1/t/ 



W^i^^(> q,T;<y,Tf) = ^erfc 



M(g) 
V2t 



l + erffm) 



\V2TfJ 



5 STATISTICAL CHARACTERISTICS OF 

DARK MATTER STRUCTURE ELEMENTS 

Results obtained above allow us to find the approximate 
PDF and other characteristics of structure elements. The 
structure element with a size (mass) m is defined as a pan- 
cake with the size m formed at a moment t that not merged 
with any other pancake. 
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5.1 Merging of dark matter structure elements 

As before the characteristics of structure element at a mo- 
ment T are expressed through characteristics of initial per- 
turbations. The approximate expression for the PDF of 
structure element can be written in a form similar to the 
known equation of coagulation: 



poo rn 

Jo Jo 



dxNcr{x, T)Nc{m—x, x, y, t 



(5.1) 



dWmerg{x,1 



,y) 



-Ncr{m, t) 



dx 



dyNc{x,m,y,T) 



dy 



dWmerg{x,m,y) 

dy 



The functions Ncr,Nc & Wmerg are given by (3.6), (3.15), 
(3.17) & (3.18). Here the first term describes the formation 
of two pancakes with sizes x k, m — x and a separation y 
and their merging to a pancake with the size m while the 
second term describes merging of the pancake of size m with 
another pancake. If the mass exchanged during merging is 
incorporated into the forming pancake then the first term in 
(5.1) has to be appropriately changed. 

Here, as the first step of investigation, we will use the 
simpler approximate approach based on the survival proba- 
bility of a pancake with size m to avoid merging with larger 
pancakes with sizes x> m. For the more interesting case of 
smaller pancakes with qo ^ m < x <1, the most prob- 
able process is the formation of two pancakes with sizes 
m & X > m at a small separation |-Dsepl < 0.5(a; + rn) that 
means, as follows from (3.1), formation of one structure cle- 
ment with a size x. This process is described by the second 
term in (5.1). Because in this case the probability of merging 
quickly decreases for larger separations |-Dsep| > 0.5(rn, + .x) 
and the function Wc{rn, x, Dsep) weakly depends on Dsep, 
we will use the approximate expression for the probability 
of merging, Pmrg, 



,(m,T) « erfc 



rV2 



m < 1, 



(5.2) 



(5.3) 



(5.4) 



^^/2 

In spite of the approximate character of this approach, it 

allows us to obtain reasonable estimates of the expected effi- 
ciency of merging and of the large scale bias. As it is directly 
seen from (5.4), for fi{m) < r, 

Psurv{m,T) oc IJ.{m)T~^, 

what characterizes the impact of pancake merging. 



and for the survival probability, Psurv , 

Paurvim, r) W 1 - Pmrg{m, t), 



Psurv{m, t) w erf ( ) , m < 1. 




2.0 



Figure 1. The PDF of structure elements, Nstr(rn,T), the mass 
distribution function, mNstr{m,T), and the fraction of com- 
pressed matter, /(> m), vs. the masses (sizes) of structure el- 
ements, m/mm, are plotted for five moments of time: t = 0.1 

(solid line), r = 0.3 (dashed line), t = 0.5 (dot-dashed line), T = 
0.7 (dot-dot-dot-dashed line), r = 2 (long dashed line). 



5.2 Statistical characteristics of structure 
elements 

With this survival probability the approximate PDF for the 
structure elements, Netr{q,T), can be written as follows: 



Nstr{m,T) oc Psurv {m,T)Ncr{m,T). 



(5.5) 



These relations allow us to obtain also the approximate mass 
distribution function, N^^\m,T), characterizing the distri- 
bution of compressed matter over the structure elements. 
For the more interesting case m <Cl we have 



Nstr{m,T) 



24 



17V27rT dm \ V2r / V 



< m >w 8r 



(5.6) 



In the general case this function can be obtained rmnicri- 
cally. The mean cumulative comoving linear number density 
of structure elements, nstr{> tUjT), is 



nstr{> m,T) 



, dn(> m,T) ^ , , _ „, 

dm ; Psurv(m,T) (5.7) 



dm 



n{> m,T)Psurv{'rn,T) + / dmn{> m,T) 



,dPs 

dm. 



Some of these functions are plotted in Figs. 1, 2 & 
3. Fig. 1 shows the PDF, Nstr(m,T), the mass function, 
mNstr(m,T), and the currmlative mass function, /(> m), 
for structure elements versus the m/mm for five time mo- 
ments. Here mm is the median mass of structure elements 
defined by the condition /(> mm) = 0.5. At small r the 
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1.0 



0.1 ^ 




0.0 0.5 



1.0 1.5 2.0 2.5 



10.0 ^ 




0.4 0.6 

l-f(>m) 

Figure 2. Top panel: the mean comoving linear density of a 
pancake, lonstr(> ™, t), vs. masses m/mm are plotted for five 
moments of time: r = 0.1 (solid line), r = 0.3 (dashed line), t 
= 0.5 (dot-dashed line line), r = 0.7 (dot-dot-dot-dashed line), 
T = 2 (long dashed line). Bottom panel: the mean separations of 
pancakes, Dstr{> m,T)/lo, vs. the fraction of compressed matter 
/(> m), for the same time moments. 



distribution function is rapidly changing duo to the merging 
of smaller pancakes. Later the evolution is slower and it is 
sustained by growth of median mass due to the progressive 
matter concentration within richer structure elements. 

The same effect is clearly seen in Fig. 2 where the mean 
linear number density of structure elements, lonstr{> rn,T), 
is plotted for the same time moments. Rapid merging of 
low mass pancakes and formation of more massive pancakes 
changes the shape of this function at m <C mm and m 3> 
rrim while at m --^ rrim evolution is very slow. The function 
DstriJ) plotted in Fig. 2 (bottom panel) is similar to that 
found in simulations (DFGMM). 

Time dependence of the mean, < m >, and median, 
mm, masses of structure elements is plotted in Fig. 3 (top 
panel) together with fits 



< m(T) >« 9tV(1 +4t), r<1.5, 



mm(r) w 0.03 + 2.2r, t < 1.5. 



(5.8a) 



(5.86) 



Faster growth of < m(r) > in comparison to (3.7) is caused 
by merging of pancakes as it is described by (5.5). In Fig. 3 
(top panel) the Dstr{mm,T) is plotted together with the fit 



Dstr{mm,T) ^ 1.5 • 



T + 0.07 



0.0 




1.0 



0.1 



1.0 



(5.8c) 



Figure 3. Top panel: Time dependence of the mean mass, < m >, 
(solid line), median mass, irim, (dashed line) and the mean separa- 
tion Dsep(iyim) /Iq- The fits (5.8) are plotted by dotted lines. Bot- 
tom panel: Evolution of the mean comoving linear number den- 
sity of structure elements, Ustrlo, vs. r for four threshold masses 
m/lo > 0.05 (solid line), mjlo > 0.1 (dashed line), m/lo > 0.3 
(dot-dashed line), m/lo > 0-5 (dot-dot-dot-dashed line). A power 
law r"" '* is plotted by a dotted line. 



Bottom panel in Fig. 3 shows evolution of the lin- 
ear number density of structure elements, Ion, caused by 
their exponentially fast formation at smaller r, and later, 
by successive merging which is well fitted by the expres- 
sion Dstr{> g, t) oc t~° * what is slower then that found in 
DFGMM. 



5.3 Parameters of structure in observed 
catalogues 

In previous Sections the formation and evolution of struc- 
ture was described in a comoving space. But in observed 
catalogues the redshift position of galaxies along the line 
of sight is used. This difference distorts the parameters of 
observed structure with respect to the theoretical expecta- 
tions (see, e.g., Melott et al. 1998; Hui, Kofman & Shandarin 
1999). To take into account this distortions the theoretical 
relations given above need to be modified. 

In observational catalogues the distance to a galaxy is 
defined by its observed velocity which can be found from 
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(2.1) 



Hjz) 
1 + z 



[Qi - p{z)Biz)Siici)], 



(5.9) 



0- 



1 + z dB{z) 



B 



dz 



1 </?<2, 



where H{z) is the Hubble constant. Therefore, for observed 
catalogues all relations obtained above will be valid after 
replacement 



t{z) — > Tv{z) = t{z){I3{z) cos ip + sini^), 



(5.10) 



where ip Is a. random angle between the hne of sight and 
the directi on of pa ncake compression. It is evident that 1 < 
Tv/t < yjfp + l, and for p = 7r/2 — r. This means 
that the observed structure parameters will be randomly 
increased with respect to those found in the comoving space. 

The value of this distortion depends on the eosmological 
parameters and redshift and is superposed with distortions 
generated by the structure disruption, selection effects, and 
other random factors. The well known example of such dis- 
tortion is the effect of so called "finger of God" . The direct 
comparison of simulated structure parameters in comoving 
and redshift spaces performed in DMRT demonstrates a 
moderate dependence of the main structure properties on 
these factors. But some properties of the structure such as 
characteristics of galaxy distribution within walls are differ- 
ent in comoving and redshift spaces. More details can be 
found in Sec. 9 and DMRT. 



6 FORMATION OF VOIDS 

The same method can be used to solve the complementary 
but more complicated problem of void formation. Accord- 
ing to the Zel'dovich theory the probability to find a 'void' 
with a size r > is identical to the probability that the 
inequality 



(1 + z)r = qi-q2- tASi2 > (1 + z)rv 



(6.1) 



holds under the conditions that pancakes have been formed 
near the boundary points qi<kq2, and in the absence of any 
pancakes between points qi&iq2- This is equivalent to the 
probability that the inequality 



A5'i2 < [qi2 -{1 + z)rv\/T 



(6.2) 



holds under the same conditions. 

More accurately we can consider the under dense re- 
gions bounded by pancakes with masses mi&7Ti2 (or with 
masses exceeding the threshold mass, m\,m2 > rrithr) in 
the absence of any pancakes with masses m > rrithr be- 
tween points gi&g2. The methods discussed above and sup- 
plemented by the void probability function technique could 
be used, in principle, to obtain such probability but the 
derivation becomes very cumbersome due to many addi- 
tional conditions. Numerical methods similar to that used 
by Sahni et al. (1994) could be more suitable for such an 
investigation. 

Let us remind however that the parameters of under 
dense regions are closely linked to the spatial distribution 
of initial gravitational potential what allows to specify con- 
veniently the large scale perturbations. This link is clearly 



demonstrated by the adhesion approach (see, e.g., Shandarin 
& Zel'dovich 1989) and was discussed in details in Paper I. It 
can also be described with the technique considered above. 
Thus, for example, the large scale modulation of the pancake 
distribution by the spatial distribution of this potential can 
be described as a modulation of the effective 'time' moment 

Teff- 

To illustrate this statement we can consider the pancake 
formation for a given potential difference, Ai^, between two 
points with comoving coordinates xi & xi + Dcm- In this 
case the distribution function of the pancake with a size 
m <C DcM at the point with a comoving coordinate x\ < 
y < xi+ Dceii is given by (3.6) & (5.5) with 

3_ ^ i ( 1 - K^(^Dcell,y) — ] , 0<y< Dcell, (6.3) 



2Go(Dc 



= Gi2{y) - Gi2{Dc 



Here is the dispersion of gravitational potential difference 
at the points with the separation Dceii- Functions Go&iGi2 
are introduced in Appendix I. 

As is seen from (6.3) the efficiency of pancake formation 
depends directly on the coordinate y and the potential dif- 
ference, A4>{Dcen). It is amplified in regions with Ktt,Acj) <0, 
and depressed in regions with it^Atj) >0. 

The same large scale modulation modifies the Euler po- 
sition of pancake and instead of Eq. (3.16) we obtain after 
averaging of A<p{m)/m (under the condition m <C Dceii) 

A(j> 



<{l + z)rpan >=y + rGriy, Dcell)- 



(T4> 



(6.4) 



Gr{y,D) = yGUy) + {D- y)Gi2{D - y), 
and for the mean Eulerian linear number density of pancakes 



rieiv) = 



dy 



1 +TGn{y, Dcell) 



(6.5) 



G„(2/,D) = Gi2{y)-G^2{D-y)+y''G2z{y)-{D-yfG2z{D-y). 

These expressions describe the matter infall into wells 
of initial gravitational potential and clearly demonstrate the 
close comiection of under dense regions with the spatial dis- 
tribution of initial gravitational potential. More detailed dis- 
cussion of this effect can be found in Paper I. 



7 LARGE SCALE BIAS AND 

CHARACTERISTICS OF BIASED COLD 
MATTER DISTRIBUTION 

7.1 Large scale bias 

In observations we have to deal with the galaxy distribution 
which is, probably, biased relative to the spatial distribu- 
tion of dark matter in the considered range of scales. The 
qualitative physical ground of such large scale bias is trans- 
parent and was discussed already ten years ago by Dekel 
& Silk (1986), and Dekel & Rees (1987). It can be caused 
by the joint action of large scale perturbations and reheat- 
ing of the baryonic component of matter. The interaction 
of large and small scale perturbations shows up as an ex- 
cess of low mass pancakes and acceleration of their growth 
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within richer waU-hke elements and is seen in simulations as 
a strong concentration of high density peaks within filaments 
and richer walls. It results in an excess of compressed cold 
gas inside "proto" walls before reheating. But after reheat- 
ing further compression and cooling of gas are inhibited duo 
to the growth of entropy of uncompressed gas. This excess of 
cold, low entropy baryonic matter, reached at the reheating 
redshift, z = Zh, can be considered as the excess of galax- 
ies within richer walls that is the large scale bias in spatial 
galaxy distribution. 

In Paper I the modulation of the rate of pancake forma- 
tion by the large scale perturbations of initial gravitational 
potential was discussed. Here we consider the direct inter- 
action of earlier and later pancakes, that allows us to obtain 
more reliable estimates of the large scale bias. With the tech- 
nique developed above we can consider only the direct inter- 
action of two population of pancakes, namely, those created 
at z = Zh, and at z = 0, and we have to neglect the mutual 
influence of intermediate population of pancakes accumu- 
lated by the walls. Such an influence can in principle be 
considered but the description becomes very cumbersome. 
So, our estimate presented below is actually the lower limit 
of the possible large scale bias. Perhaps, this problem can 
be studied in simulations which take into account all impor- 
tant factors together (see, e.g., Sahni et al. 1994; Cole et al. 
1998). 

For the warm dark matter model discussed in Paper I 
the mass of DM particles restricts the minimal size of created 
pancakes to the correlation scale of initial density pertur- 
bations, Tc. For the CDM power spectrum considered here 
there are no natural limits and, in principle, pancakes with 
arbitrarily low mass can form. Formally, it is due to the di- 
vergence of higher moments of power spectrum and, because 
of this, the typical scale rc is zero. In this case, however, 
the typical scale <7o introduced by (3.4) plays similar role 
and discriminates the process of earlier low mass pancake 
formation, described through the deformation tensor, and 
later formation of structure elements discussed above. The 
typical mass of such low mass pancakes given by the expres- 
sion (3.4) exceeds the estimates of minimal baryonic mass 
formed before the reheating, Mmin ~ W^Mq (Tegmark et 
al. 1997), and the value go can be accepted as a natural limit 
of pancake size. Further evolutionary history of low mass 
pancakes is uncertain, and this problem should be consid- 
ered separately. 

Thus, for the CDM power spectrum we assume that: 

(i) All baryons accumulated by structure elements with a 
size m > go, at the redshift z — z^, t — t^, arc incorporated 
into observed 'galaxies'. 

(ii) The formation of such clouds is interrupted at the 
redshift z = Zh, t = th hy the instantaneous reheating of 
the uncompressed gas. 

Both assumptions are very restrictive but they can be used 
for rough estimates of the efficiency of this mechanism of 
generation of the large scale bias. The moment of reheating 
is, probably, restricted to 7 < zu < 15. The lower limit is 
imposed by the observational constrains. The upper limit 
is imposed by the very fast growth of the inverse Compton 
cooling for larger redshifts. 

This means that for the models with fim ~ 0.3 - 1 the 
reheating probably occurred at Th = T{zh) ~ (0.15 - 0.3)to. 



7.2 Statistical characteristics of biased cold 
matter distribution 

The biasing factor on the LSS and SLSS scales can be es- 
timated in the same manner as the probability of merging 
(5.2). Now we are interested in the mass of cold matter col- 
lapsed at T = Th which was accumulated at the moment r 
by a pancake with a size D S> nicid- The approximate PDF 
for this mass can be written as 



Ncid{mcid, Th) w 2iVc ( rricid, Th; D, r. 



Ncr{mcld,Tl), 



(7.1) 



_1_ 



1 

Th 



Kl{D) 



{2 + D){l + D)' 



These relations show that, as it was discussed in Paper I, 
the formation of earlier pancakes within larger structure el- 
ements is accelerated, what is seen in (7.1) as a shift of the 
effective evolutionary 'time', ri, relatively to the real re- 
heating 'time' moment, Th- The corresponding mean 'mass' 
of such pancakes can be found as 



< mcid{Th) >w 4ri w 4r^(l -|- 1.5r^/r), 



D<1. (7.2) 



The dimensionless biasing factor 
< rricoid > 



bbias 



<m> 



- 1 



1.5rh/r, 



(7.3) 



is about of 1/5 - 2/5 and weakly depends on the size D of the 
considered structure element. The relation (7.1) describes 
also the excess of larger pancakes induced by the influence 
of larger structure elements 



Ncid{mcid,Th) 

Ncr{fJ',Th) 



Th 



(7.4) 



As is seen from (7.1) the bias decreases, for t ^ Th- This 
means that the indirect multi step interaction becomes spe- 
cially important and the intermediate population of pan- 
cakes formed at < r* < r with a size x* < D and 
accumulated by a structure element with a size D, at the 
moment r, will enhance the bias and increase the concen- 
tration of cold matter within larger pancakes. 

These results are consistent with the conclusions of Pa- 
per I, Cole et al. (1998) and Sahni et al. (1994) and demon- 
strate the sensitivity of the considered bias to the detailed 
characteristics of the process of galaxy formation. They show 
that the over and under dense regions found now in the ob- 
served galaxy distribution can, probably, be related even at 
high redshifts to a spatial modulation of sizes of the earlier 
pancakes. The larger pancakes are more closely associated 
with over dense regions whereas the under dense regions are 
occupied preferentially by low mass pancakes. 

Summarizing, we can infer that the influence of large 
scale bias generated by a combined action of reheating and 
large scale perturbations can be important and can distort 
characteristics of the observed galaxy distribution relative to 
the same characteristics of DM distribution. Under suitable 
conditions the fraction of baryons accumulated by galaxies 
can be suppressed in low massive pancakes and under dense 
regions and increased in richer structure elements. 



© 0000 RAS, MNRAS 000, 000-000 



12 Demianski & Doroshkevich 



8 DYNAMICAL CHARACTERISTICS OF THE 
PANCAKES 

In this section we will find two dynamical characteristics of 
pancakes, namely, the velocity of pancakes as a whole and 
the dispersion of matter velocity within pancakes caused by 
compression of pancakes. 

The velocity of an infalling particle, Vi, can be found 
from (2.1) as 



1 + z' 



0{z)B{z)Si{cD], (8.1) 



0=1- 



1 + z dB{z) 



B dz ' 

where H{z) is the Hubble constant. 

The functions /3{z) can be found from (2.2) and (2.3) 

as 

1 f^m 



B^z) 



1 + i.2a„ 



B{z) 



1 

1 + i.5n„ 



+ = 1, (8.2) 



< 1, = 0. (8.3) 



For fim = 1 both expressions give /? = 2. 

Using the general relation for the velocity of a particle 
(8.1) we can find the velocity of a pancake as a whole through 
the momentum conservation law, we have 



1 

qi2 



v{q)dq ■■ 



Hjz) 
1 + z 



Qcent - I3B{Z 



A<^12 



qi2 



A012 = </>(gi) - ^(92), gcent = (91 + 92)72. 



(8.4) 



The pancake position, Vpa-, 
liar velocity of pancake is 



was given by (3.16). The pecu- 



v,ar.-H{z)r,.^ = ^B{(i-l)^'^'' 



^, ■ (8.5) 
l + a; ' qii 

The mean peculiar velocity, < u >= 0, and its dispersion. 



2, , -ff^(^) rr,, N ii2d2/ n2Go(912) 

I12 



The distribution function of peculiar velocity is Gaussian as 
before. In dimensional variables, aX z = 0, and for 912 <S 1, 
we have: 

o-„(0) = Ho[l3{0) - l]as/V3 = Ho[f3{0) - 1]/oto. (8.7) 

These relations show that for smaller 912 the pancakes ve- 
locity is close to the mean velocity dispersion and decrease 

X /2 

oc qj^2 ) for qi2 > 1. The vortcxlcss character of the initial 
velocity perturbations is clearly seen in eq. (8.4), (8.5). 

Variations of the pancake velocity with the transversal 
coordinates can distort the shape of pancake and, in princi- 
ple, can destroy the pancake. The impact of this factor can 
be characterized through the expression (8.4) as a correla- 
tion of pancake velocity at different transversal coordinates. 
This correlation is described by the normalized function r„ 
(see Appendix IV), and for 912,934 <C 1, we have 

rv = qi2 • q34/9i2934, (8.8) 

what demonstrates the regular character of velocity varia- 
tions. For larger 912 and 934 the velocity dispersion (8.6) 
decreases and therefore this factor becomes less important. 
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Figure 4. The velocity dispersion, (T„, vs. the size of "proto 
walls", q = R-wall, for the redshift z = 0, for five cosmological 
models with h = 0.5, Qm = Ij^a = (solid line), and h = 0.7, 
Clm = 0.5, Ha = 0.5 (dashed line), dm = 0.5, Ha = (dot- 
dashed line), Qrn = 0.35, Qa = 0.65 (dot-dot-dot-dashcd line), 
f^m = 0.3, £1a = (long-dashed Unc). Dotted lines show observa- 



tional estimates of the size of " proto walls" 
and the velocity dispersion it„ ~ 400km/s. 



ill 



' 30h-iMpc, 



The velocity dispersion within a pancake with a given 

mass q = qi — q2, caused by the matter infall into pancake, 
is defined by the energy conservation law. The velocity of a 
particle with the coordinate 93 relative to the pancake as a 
whole is 



i{q3,z) = v{q3,z) - u{z) 



1.9) 



Hjz) 
1 + z 



93 - qcent - l3{z)B{z) S{q3) 



A(?!)i2 
912 



Using the conditional mean values and dispersions of ran- 
dom functions 5(93), and 1^(91) — 4>{q2), listed in Appendix 
I, we obtain the required velocity dispersion 



(1 + ^)' 2 
-cr,. 



H^{z) 



12 



mz)-lf+2P{z)h{q)- 



f{z)f2{q) + T'l3\z)f3{q)]. 



(8.10) 



The functions /i, /2, & /s are given in Appendix IV. The 
velocity dispersion (Jv{q) is plotted in Fig. 4 for a few cos- 
mological models, for z = 0, h = 0.7. 

For 912 < 1 the functions /i,/2&:/3 are small, /i < 
0.06, /2 < 0.13, and fz < 0.2, and the contribution of the 
three last terms in (8.9) becomes important only when /? ~ 1 
that is for cosmological models with a lower matter density 
at redshift z < 1. For poorer pancakes when the first term 
dominates, the simple approximate relation 



ay^H{z){l + zy\f3- 



(8.11) 



can be used for estimates of the oxpoctod velocity dispersion. 
For such pancakes the distribution function of av is similar 
to (4.5). 
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9 COMPARISON WITH SIMULATED AND 
OBSERVED STRUCTURE PARAMETERS 

Some characteristics of large scale galaxy distribution can be 
extracted from available catalogues of galaxies and clusters 
of galaxies and catalogues of absorbers in spectra of quasars. 
Some characteristics of large scale DM distribution can bo 
derived from available simulations. All these characteristics 
are distorted by selection effects, bias between spatial distri- 
bution of DM and galaixies, and other factors what compli- 
cates the direct comparison of characteristics obtained with 
different methods and data bases. 

Nonetheless, the comparison of even so distorted char- 
acteristics with the expected characteristics discussed above 
can be interesting as it allows one to compare various esti- 
mates and illustrates their sensitivity to basic cosmological 
parameters and other factors. More detailed discussion can 
be found in DMRT & LCRS2. 



9.1 Observed and simulated peirameters of the 
structure 

Analysis of galaxy distribution in the Las Campanas Red- 
shift Survey (Shectman et al. 1996) shows that the richer 
and poorer structure elements can be assigned to wall-like 
and filamentary populations, but the accurate demarcation 
of these subpopulations is problematic. Both subpopulations 
accumulate ^ 40 - 60% of galaxies. Filaments fill the gaps 
between the walls and form a random network with the mean 
cell size Df' ~ 10 - 12^-^ Mpc at « = 0. The richer walls 
can be formed due to an anisotropic compression of matter 
within slices with a typical size of R'^aii ~ 20 - 25/i^^Mpc 
and with a typical separation of D^^n ~ 50 - 60/i.~^Mpc 
(LCRSl & LCRS2). The walls are significantly disrupted 
by the small scale clustering of galaxies (see, e.g.. Fig. 5 in 
Ramella et al. 1992). The velocity dispersion within such 
wall- like elements was estimated as ai°''^^ ~ 350 - 400 
km/s (Oort 1983a). The bulk velocities of galaxies are now 
estimated as a"'" ~ 400km/s (see, e.g., Dekel 1997). The 
observed surface density of structure elements is heavily dis- 
torted by the selection effects and small scale clustering and 
is equally well fitted to a power law, exponential distribu- 
tion, and functions discussed in Sees. 4 & 5. 

The structure in the DM distribution was simulated and 
analyzed by DMRT for the SCDM model with 0^=1, h= 
0.5, ACDM model with 0^=0.35, nA=0.65, h = 0.7 and 
OCDM model with $1^=0.5, h = 0.6. For all these models 
the normalization to the two year COBE data were used. 
These simulations reproduce the structure with the main pa- 
rameters similar to those found in observations. The analysis 
shows that walls accumulates ~ 40% of DM with a mean size 
of "proto walls" R'^^ 15 - 20/i"^Mpc and with a mean 
separation of walls and filaments D^J^i ~ 50 - 70/i~^Mpc, 
Df"" ~9 - ll/i-^Mpc, respectively. 

In all the models the simulated velocity dispersion 
within walls, uj*'", is the same along all three principal di- 
rections of walls and it is generated by the wall curvature 
and disruption into a system of high density clumps. The 
walls are seen usually as ordered sets of irregular high den- 
sity clumps connected by lower density bridges. The degree 
of wall disruption depends on the code used and reached 
resolution, and it is more conspicuous in the models with 



larger Qm- In the redshift space this inner structure is partly 
eroded and characteristics of DM walls are similar to that 
found for the observed galaxy distribution. 

The isotropy of simulated velocity dispersion (t„ con- 
firms the essential influence of small scale clustering on the 
properties of both observed and simulated structure ele- 
ments. The instability of a thin compressed layer was an- 
alyzed in the linear approximation by Doroshkevich (1980), 
and Vishniac (1983), and was recently simulated by Valinia 
et al. (1997). The clustering rate and parameters of formed 
clusters depend on the density of compressed matter and 
properties of transversal velocity field. 

9.2 Comparison of simulated and expected 
structure properties 

Theoretical parameters of DM wall-like structure elements 
are listed in Table I for the median mass m = nim and for 
Ic/lo =0.056. The variations of Rwaii listed in Table 1 are 
moderate and, in the range of precision of our approximate 
consideration, these values axe similar to those observed and 
simulated. The differences become more significant only for 
the OCDM models. The value of Dstr is sensitive to the as- 
sumed mass threshold (see Figs. 2 & 6) and can be adjusted. 

The simulated structure elements are also distorted by 
the small scale clustering but owing to the rich statistics of 
such elements a more detailed quantitative comparison of 
simulated and expected properties of the structure can be 
performed. To do this we will consider the PDFs (4.4) and 
(4.5) for the surface density of filaments and pancakes, the 
PDFs for velocity u and velocity dispersion it„ and direct 
estimates of important parameters 



> 



HoloW - 1) ' 



mv/lo 



2V3 < g-g > 
HoloiP - 1) ' 



(9.1) 



following from relations (8.7) and (8.11). 

To characterize the filamentary and wall-like subpopu- 
lations of structure separately the analysis was performed 
both in the comoving and redshift spaces for three subsam- 
ples of structure elements. The first contained all particles, 
the second and third incorporated the richer and poorer 
structure elements with Nmem > 200 and Nmem < 200, 
respectively. Here Nmem is the number of DM particles 
within a structure element bounded by the threshold den- 
sity nthr/ < n >=1.7, 1.1 & 1.1 for the SCDM, OCDM 
and ACDM models. In all the cases the second subsamples 
accumulates ~40 - 45% of particles. 

To obtain the required PDF a set of rectangular sam- 
pling cores was prepared and the number of DM particles 
in the intersection of separate structure elements with these 
cores was taken as a characteristic of surface density of struc- 
ture elements. For each of the cluster its velocity u and veloc- 
ity dispersion cr„ along the core were also found. To depress 
the impact of small scale clustering the analysis was per- 
formed for the core size Lcore = 10/i~^Mpc, for the second 
subsample, and for Lcore = 4/!,~^Mpc, for the first and third 
subsamples when the core size is restricted by the separation 
of filaments. The random intersection of cores and structure 
elements generates significant excess of poorer clusters. To 
depress this effect poorer clusters were rejected and the trun- 
cated PDFs were considered. Even so, the number of clusters 
used was ~ 3 500, for the second, and ~ 10 000, for the first 
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Table 1. Main expected parameters of the SLSS. 
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Figure 5. The simulated PDFs of surface density of structure ele- 
ments, Nstri'm), vs. m/ < m > for the full sample (top panel) and 
subsamples of richer (middle panel) and poorer (bottom panel) 
structure elements for the ACDM model. The fits (9.3) and (9.4) 
arc plotted by solid lines, the power fits are plotted by dashed 
lines. 
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Figure 6. The simulated PDFs of velocity dispersion of struc- 
ture elements, Nstr{<^v), vs. ay/ < (Jv > for the full sample (top 
panel) and subsamples of richer (middle panel) and poorer (bot- 
tom panel) structure elements for the ACDM model. The fits (9.3) 
are plotted by solid lines. 



and third subsamples. The random orientation of structure 
elements and cores can increase the measured surface den- 
sity (up to two times for homogeneous matter distribution, 
Kendall & Moran 1963), and reduces the measured velocity 
of structure elements along the core, Mc, by a factor of \/3 
and we have 

Wc = MCOS(^, T„ = \/3t„c. (9.2) 

To compare the simulated and theoretical PDFs the two 
parameters fits were used. For the second subsample the 
distribution of surface density was fitted to the expression 

N^ = ^e-^-er!(V^), = (9.3) 

which reproduces approximately the theoretical relations 
(4.5) with the correction to the pancake merging discussed 
in Sec. 5. For the first and third subsamples the expression 

Nf = ^e-f[l + erfi^)], = ^ (9.4) 



reproduces approximately the theoretical relations (4.4) and 
gives better fits. Here parameters by, &Lbf describe the trun- 
cated character of fitted PDF whereas parameters ayj &i a/ 
provide its normalization. For velocity of clusters, u, the 
PDF is found to be well fitted to the expected Gaussian dis- 
tribution, the expression (9.3) with parameters o„ & b„ fits 
the velocity dispersion cr„. 

For ACDM model these PDFs are plotted in Figs. 5 
and 6 together with the best fits. It turns out that different 
PDFs verify the accepted discrimination of two populations 
of structure elements. For the second subsample the theo- 
retical estimate of mean wall size 

< Trip > /lo w 8(0.5 + l/7r)ro » 6.55ro , (9.5) 

links the parameters To and < rUp >. For the ACDM and 
OCDM models the parameters Tu (for all samples) and < 
mv/lo >, Tv = \J < ruv/lo > /6.55, < rup/lo > <k Tp = 
y'<~m^7^o3~76^ (for second subsample) listed in Table 
I are consistent with the simulated amplitude. For the first 
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and third subsamples the distributions plotted in Fig. 6 and 
corresponding parameter < m„/Zo > characterize the typical 
mass of dominant numerous poorer elements and it is about 
half of that for the second subsample. 

The formal precision of these estimates can be taken 
as ~ 7 - 10% (precision of fits) but a real precision can be 
estimated by the comparison with the value of the input pa- 
rameter To also listed in Table 1. Differences between results 
obtained in comoving and redshift spaces do not exceed the 
reached precision. For the SCDM model the expected and 
reconstructed parameters differ by a factor of ~ 1.5 what 
can be partly caused by the strong wall disruption. More- 
over, this simple description is correct for m/lo <1 whoroas 
more cumbersome general relations describe the properties 
of richer walls which dominate the SCDM model. 

In all these cases the simulated mass distribution can 
be equally well fitted to a power law with the exponent k ~ 
1.5-2. Such distribution is similar to that described analyt- 
ically for the matter concentration within a set of clusters 
with a surface density acis oc r~'^ (see also discussion in 
Miralda-Escude et al. 1996). In this case the measured PDF 
and mass function are also expressed by the power law 

dWcis oc adsrdr oc a~il^'^dacis- (9-6) 

Of course, profiles of separate irregular clusters can vary 

over a wide range and such description and interpretation 
illustrate only the important role of small scale clustering. 



9.3 Evolution of DM structure 

Results obtained in Sees. 3, 4 & 5 show that during the 
most interesting period t ^ Ic/lo evolution of DM structure 
shows features of self-similarity and important characteris- 
tics of structure can be expressed as functions of the pa- 
rameter /i(g)/r. This approximate self-similarity is caused 
by the Zel'dovich approximation and occurs for any dis- 
tribution function of initial perturbations and any power 
spectrum. It becomes more transparent for simple structure 
functions, at go g <1 (Appendix I), typical for the CDM 
like initial power spectrum, but becomes more cumbersome 
for T ~ Ic/lo <^ 1, and t > 1, when the influence of the scales 
qo&iLo, introduced by (1.6), (1.7) becomes important. This 
self-similarity allows us to characterize the structure evolu- 
tion by the time dependence of several typical parameters 
such as the mean and median masses, < rn > & rn,n- 

This approximate self-similarity is partly violated due 
to the evolution of pancakes, discussed in Sec. 4, and small 
scale pancake disruption as they are not described by the 
Zel'dovich theory. Nonetheless, the observations and simu- 
lations show that in a universe dominated by cold DM parti- 
cles general properties of matter distribution are quite simi- 
lar at high and small redshifts (see more detailed discussion 
in Govornato ct al. 1998; Jenkins et al. 1998). Self-similarity 
of structure evolution was previously discussed for the scale- 
free power spectra (sec, e.g., Efstathiou et al. 1988). 

The evolution of observed structure formed by galaxies 
can however be far from the self-similar evolution of DM 
structure due to the small and large scale bias. Nonetheless, 
as it was discussed above, the main parameters of DM struc- 
ture can be compared with observed galaxy distribution. 
Thus, in particular, the observed separation of filaments, at 



a = 0, is about 4-6 times smaller then the typical separa- 
tion of walls (see, e.g., Efstathiou et al. 1988). Such result 
can be attributed to pancakes formed by compression of a 
slice with thickness m/ ~ (0.1 - 0.05)mm « 1 - 2/i~^Mpc 
at redshift z ^ 3 5. In smaller pancakes the formation of 
gala^xies can be suppressed and they can be associated with 
weak Ly-a absorbers observed far from galaxies (Morris et 
al. 1993; Stocke et al. 1995; Shull et al. 1996). 



10 SUMMARY AND DISCUSSION 

In this paper we continue investigations and statistical de- 
scription of the process of structure formation and evolu- 
tion initiated in our previous papers (Buryak et al. 1992; 
Demianski & Doroshkevich 1992, 1997; DFGMM; Paper I) 
and based on the nonlinear theory of gravitational insta- 
bility (Zel'dovich 1970). The new elements discussed above 
are the approximate statistical description of the nonlinear 
structure evolution manifesting itself as a successive mat- 
ter concentration into more and more massive structure ele- 
ments and the estimates of the large scale modulation of the 
spatial distribution of luminous matter relative to DM and 
baryonic components. We show that, as has been discussed 
in Sec. 9.4, the evolution of DM structure demonstrates some 
features of self-similarity and the main characteristics of the 
DM structure can be expressed through the structure func- 
tions of initial perturbations and through parameters of the 
initial power spectrum and cosmological models. It is shown 
- in accordance with simulations - that in a low density 
models the nonlinear evolution occurs on the scale --^ 20 - 
25/i~^Mpc and results in formation of wall-like component 
of the structure of the universe. 

Results discussed in Section 3 show that, for the CDM- 
like power spectra, the high nonlinear matter compression 
into low mass pancakes and the percolation take place al- 
ready at T « Ic/lo <d. Later evolution leads to a rapid 
growth of typical mass of DM pancakes. The theoretical 
description discussed above can be directly applied to the 
interpretation of deep pencil beam galactic surveys and ab- 
sorption spectra of quasars, what allows us to consider to- 
gether rich observational data accumulated by these meth- 
ods for both small and large redshifts. For larger 3D catalogs 
and simulations this description can be applied to results 
obtained with the core-sampling analysis as it was demon- 
strated in Sec. 9.2 . 

Comparison with simulations shows that the velocity 
dispersion of matter compressed within the wall-like ele- 
ments is more sensitive to the cosmological parameters, large 
scale modulation of spatial galaxy distribution and proper- 
ties of DM component. The small scale clustering generates 
the essential difference between the expected and simulated 
velocity dispersion and allows us to discriminate the cosmo- 
logical models with respect to the matter density, Qm, and 
the composition of DM. This clustering is smaller for the 
models with Qm < 0.5, but for models with fim < 0.3 it is 
more difficult to reproduce the sizes, separations of observed 
walls, and the high matter concentration within the SLSS. 
The observed structure parameters arc best reproduced for 
the models with ^Imh « 0.2 - 0.3. These values are consis- 
tent with estimates obtained both from the simulations of 
cluster of galaxies (Bahcall & Fan 1998; Cole et al. 1997, 
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1998) and from the observations of high redshift supernovae 
(Perlemutter et al. 1998). 

Two factors can extend the set of acceptable cosmologi- 
cal parameters. The first is a more complicated composition 
of DM component, what means that an essential fraction of 
DM can be associated with relatively hot and/or low mass 
particles (see, e.g., Colombi et al. 1996; Brustein & Hadad 
1998). The simple estimates based on the Trimain - Gunn 
relation (Trimain & Gunn 1979) show that if the recently 
formed wall-like SLSS is sensitive to the influence of low 
mass relic particles, with Mp < 3 - 5 eV, then the structure 
properties at higher redshifts are sensitive to more massive 
particles as well. This means that for the MDM models with 
more complicated DM composition the agreement between 
observed and simulated properties of the LSS and SLSS can 
be achieved in a broader range of Qm and h and, moreover, 
the analysis of observed structure evolution can specify the 
DM composition. The models including some fraction of un- 
stable dark matter particles (Turner et al. 1984; Doroshke- 
vich et al. 1989) seem also to be perspective. Unstable parti- 
cles can produce reheating by photo decay whereas the hot 
products of decay will delay both formation and disruption 
of the walls. 

The second factor is the large scale modulation of spa- 
tial galaxy distribution. The formation of LSS and SLSS 
is always accompanied by bias on the same scales, because 
both processes are caused by the same large scale perturba- 
tions and, so, are strongly correlated in space. The efficiency 
of such bias depends on many factors, such as the redshift 
of reheating period, and the formation and evolution of cold 
clouds. The rough estimates of the bias presented in Sec. 
7 show that the bias factor could be quite high (see also 
Demianski and Doroshkevich 1997; Paper I) and together 
with other factors it can essentially distort the spatial dis- 
tribution of galaxies with respect to the DM distribution. 
These estimates could be enhanced by taking into account 
the mutual interaction of pancakes with various sizes and 
moments of formation. More representative estimates of the 
bias can, probably, be obtained from simulations such as, 
for example, Sahni et al. (1994) and Cole et al. (1998). 

This mechanism of bias generation implies early, for 
z > 5, formation of the main fraction of low entropy gaseous 
clouds, that can be identified with 'proto galaxies'. The re- 
heating does not prevent further formation of DM pancakes, 
which can be identified with a population of gas clouds re- 
sponsible for weaker absorption lines observed at high red- 
shifts (Miralda-Escude et al. 1996; Hernquist et al. 1996) but 
it inhibits the formation of high density, low entropy gaseous 
clouds and delays the formation of galaxies in extended un- 
der dense regions. Further transformation of formed earlier 
'proto galaxies' into observed galaxies is a slow and compli- 
cated process continuing up to now and earlier formation of 
such clouds is not in contradiction with the observed peak of 
galaxy formation at redshifts z ~ 2 ~ 3 (Steidel et al. 1996) . 

The observed strong variations of galaxy distribution 
with respect to the DM and intergalactic gas distribution, 
Pgai/Pgas, cau be, partly, associated with this bias. Indeed, if 
in clusters of galaxies this ratio is found to be pgai / Pgas « 0.2 
(see, e.g.. White et al. 1993) then, for example, within Boots 
Void p gal I p gas —> (Weistrop et al. 1992). The existence 
of 'invisible' structure elements, which are now seen as gas 
clouds responsible for weak Ly-a absorption lines situated 



far from galaxies (« 5 - 6/i~^Mpc, Morris et al. 1993; Stocke 
et al. 1995; ShuU et al. 1996) can also be considered as an 
evidence in favor of large scale bias. 

Simulations show also existence of bias in the distribu- 
tion of DM and 'galaxies' identified with the highest peaks 
in the initial density distribution (see, e.g.. Eke et al. 1996). 
Such bias is similar, in some important respects, to that 
generated by reheating. It operates on really large scales 
(Bower et al. 1993) and can be essential for the SLSS for- 
mation as well. It results in structure composed of filaments 
and sheets in the distribution of 'galaxies' (Doroshkevich, 
Fong & Makarova 1998). Other mechanisms of bias for- 
mation discussed recently (Coles 1993; Sahni & Coles 1995; 
Tegmark & Peebles 1998), operate on significantly smaller 
scales. 

Simulations allow us to test the joint action of vari- 
ous factors on the structure parameters and, therefore, they 
now seem to offer more perspective way for detailed inves- 
tigations of the structure formation and evolution. They 
need, however, to be essentially improved in order to dis- 
criminate the spatial distribution and other parameters of 
'galaxies' and dark matter. The methods used for the mea- 
surement and description of simulated and observed struc- 
ture should be also improved as now they provide us with 
limited information about structure properties. These re- 
strictions become especially important at higher redshifts, 
where the simulated density contrast is small. Large disper- 
sions of measured values make also any comparison of sim- 
ulated and theoretical structure parameters more difficult. 
Nonetheless, such approach seems to be perspective because 
further progress in these directions can be reached. 
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Appendix I 

Correlation functions for initial perturbations. 

To investigate mutual interaction of perturbations and 
to reveal their influence on the nonlinear processes of struc- 
ture formation we can use the conventional distribution 
functions and conventional mean values. In this appendix we 
introduce a few correlation and structure functions which 
describe the relative spatial distribution of important pa- 
rameters of perturbations. 

We begin with the structure function of gravitational 
potential perturbations which characterizes correlation of 
the gravitational potential in two points qi and q2. As 
the power spectrum is a function of the absolute value of 
wave number only this structure function depends on 
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LOO 
x/l„ 



10.00 100.00 




Figure 7. Structure functions Go (top panel) and Go, Gi, G12 
& — G23 vs. Lagrangian coordinate q/lo for the SCDM power 
spectrum. 



912 = \qi — ?2| and for the perturbations of gravitational 

potential we have 

„< A012A034 > 



= Go(a;i4) - Go(a;i3) +Go(a;23) - Go(a;24), 



A012 = (?!>(qi) - 0(q2), A034 = <?i(q3)-<?i(q4). (/.I) 

Differentiation of this structure function gives other struc- 
ture functions: 



3 < A(t>{xi2)Si{X3) > 

3 < Si(xi)Sj[x2) > 



= {xi-X3)iGl{xi3)-{X2-X3)iGl{X23), 



' AijGi{qi2) + ^^Gi2(a;))|x=xi2, 



3/0 < Sk{x\)dij{x2) > 



(7.2) 



(XfcAij +XiAkj + XjAik)G2{x) + ^^^J^^'' G23{x))\x=X2i, 

where the projection operator Atj = Sij — XiXj/x^ is used 
and X = q/lo- The scale lo was introduced by (2.5). 

In the general case these functions can be expressed 
through the power spectrum, p{k), and spherical Bessel func- 
tions and we obtain 



Go(<?) 



Gi{q) 



3 rm 
3 f°° 

-.3/0 2 / P{k)y~^^^J3/2{y)dk = G'o/x, 



G2{q) = ~ 



ill r 

(27r)3/V? J, 



k^p{k)y ^'^Jz/2{y)dk = G'i/x, 



074 p°o 

^'^^^ = (27r)3/Vg X k'p{k)y-"^Jy2{y)dk = G^a^, 
where j/ = feg, a: = q/lo- 

Gi2{x) = Gi[x) + x^G2(a;) = Gq = (xGi)', 

G23(a;) = 3G2(a:) + x'^Giix) = G'12/x, (I A) 

G234(x) = 3G2(a;) -|- 6x^G3(a;) + x^G^x) = g'iV 



Gi(0) =Gi2(0) = 1, G23(0) = G234(0) = --|, (7.5) 

and the scale Ic was introduced by (2.5). 

For realistic power spectra these functions cannot be 
expressed analytically and they have to be calculated nu- 
merically. For the Harrison - Zel'dovich primordial power 
spectrum and for the CDM transfer function (BBKS), for 
X = q/lo > 6Zc/?0) 6Zc/'o 1 these functions can be 
approximated by 

Go(a:) « -!— ln(l + x + aox^) 
2ao 



1 / x\/4ao — 1 



2aoV4oo - 1 " V 2 -h a; 
Gi{x) [1 + x + ooa;^]"^, 
Gi2(a;) Gi(x)[l - aox'^], 



G23{x) w — 2G?(x)(l -|- 3ooa;^ — aoX^)/x, 



G234{x) w 6Gi(a;)(l -|- 1.3aoa;^ -|- Aa^x^ — a^x'^), 

ao = 5 (lo/Lof ^0.3, Lo«4.Uo, 
where the typical scale Lo is defined as 

sin kLo ^ 



(i.t) 



p(k)dk. 



The functions (1.3) together with fits (1.6) are plotted in Fig. 

7. The more cumbersome expression 



Gi(x) « (1 go)[l + V10 + aox^]-\ (7.7) 
G12 « G?(l + qo)-'[l - aox" + ql/^/ql+x-'] 



qo 



Giyfo[i + y/i + ml/ii < 1 



can be used for all x. In this case the functions G2(a;), G23(a;) 
and G234(a;) can be found through the relations (1.4). 

Appendix II 

Distribution function of pancEikes. 

We begin with the distribution function for the differ- 
ence of displacements at points with coordinates qi and q2 . 
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The normalization and notation introduced in Appendix I 
are used. The dimensionless 'time' r is introduced by (2.9). 

Let us consider the deformation of spherical cloud with 
a diameter q. For the case g ^ go the general deformation of 
cloud can be characterized by the 2D random scalar function 
9(6', 0) = [S(q/2) - S(-q/2)] • q/q^ instead of the deforma- 
tion tensor dik- Expansion of this function into spherical 
harmonics allows one to obtain an approximate description 
of deformation of the cloud with a reasonable accuracy. Now 
we will take into account the spherical and quadrupolc defor- 
mations only. These modes describe successfully the trans- 
formation of the cloud into filaments and sheets while higher 
order deformations relate mainly to the cloud disruption. In 
this case we can consider the deformation of the cloud in 
three principal directions, namely, x,y,&iz. The dispersion 
of displacement differences in these directions are 

crl =< [S4q/2) - S^{-q/2)f >= 2[1 - Gi2(g)], (//.I) 

2 2 2 

CTy = = 

and their correlations arc described by the coefficient rxy'- 
<Tlr.y =< [54g/2) - 54-g/2)][5,(g/2) - Sy{-q/2)] > 



g^ G2{q/V2) _ 0.25^/2 / 1 + g \ 



2 1-Gi2(g) l + q/2 \ l + q/y^J 



{II.2) 



for qo <^ q. Taking into account the approximate character 
of our consideration and because of r^y ~ r^z = Tzy <ICl we 
will neglect these correlations and will consider the deforma- 
tion along three principal axes as independent. Under this 
assumption we have for the distribution of maximal displace- 
ment difference t^S^ > A5j, > ASz and for ( = ASx/o-^ 



dC 



W2-!T 



4V27r 



erf 



V2 



C>0, (77.3) 



C<0. 



These relations show that under the above assumptions ~ 
1/8 of matter is compressed in all directions and ~ 1/8 ex- 
panded in all directions, whereas ~3/8 of matter is com- 
pressed along i&y axes and ~3/8 of matter is compressed 
along X axis and expanded along y&cz axes. This means that 
7/8 of matter is compressed at least along one of the axes. 
The influence of correlation does not distort the symmetry 
between compression and expansion and, at g < 1, increases 
the probability to find AS^; > ASy > ASz > up to 0.21. 

Let us now consider the matter compression along the 
direction of maximal compression q = qi — q2. In this case 
the dispersions of displacement differences ASi = Si{qi) — 
Si{q2) are: 



V2[l-Gi2(g)] 



This function characterizes the intersection of two particles 
before the time moment t or for the pancake separation 
> q and, in fact, {8/7)Wcr is the cumulative distribution 
function for pancakes with masses > g, for a given t. 

The standard technique (see, e.g., BBKS) can be used in 
order to find the mean comoving linear density of pancakes 
with a given ASi = q/r along a random straight line. It 
is expressed through the characteristics of the derivations of 
function Q introduced by (3.2), Qi = dQ/dqi, and for go "Cl 

<Qi>=0, cT^i =< Q? >« 6/go, (77.6) 

<^d2 =< Qi >= <^d3 =< QI 2/go 

(go is introduced by (1.7)). The required linear number den- 
sity, n(> g), is 



lon{> g) 



3<Mr> ^ f Mg) \ ^ 
16n'^aiiadicrd2ffd3 \^/2t J 



(77.7) 



VQI + QI + QI exp 



2^J2 



Qi 



dQidQ2dQ3 



where Hr characterizes the random angle of the intersection. 
For Ic/lo <C 1, and g S> go, we obtain (3.10). 

Appendix III 
Interactions of pancakes. 

In order to investigate the correlations of the pancake 
properties we need to consider the intersection of the points 
gi and g2, gi2 ~ qi — q2, at the moment ri, and points gs 
and g4, g34 = ga — g4, at the moments T2, respectively. We 
consider the simplest case when the vectors qi2 and q34 are 
situated along the same line. 

It is convenient to characterize the pancakes by their 
sizes, Di and D2, and by the separation of their central 
points, Daep. In this case we have 

qi2 = Di, g34 = D2, 
qsi = D,ep - 0.5(7)i - 732), 941 = D^ep - 0.5(7)i + D2), 

932 = Dsep + 0.5(7)i + D2), 942 = D^ep + 0.5(7)i - 7)2). 

The dispersions and the correlation of displacements 
AS{Di) and AS'(732) can be written as 



o-?(7)i) =cTii(73i), ffi(732) = fTii(732), 



(777.1) 



cT?i =< ASiASi >= 2[1 - Gi2(9)], 



(77.4) 



iTio-2 rs = Gi2(g3i) - Gi2(g32) + Gi2(g42) - Gi2(g4i), 



ai2 =< AS2AS2 >= ah =< AS3AS3 >= 2[1 - Gi(g)], 

< A5iA52 >=< AS1AS3 >=< A52A53 >= 0. 

For the distribution function (II. 2) the probability to 
have ASi > g/r is 



Wcr{q) 



1 + erf 



\V27 



(77.5) 



and rs is a symmetric function of Dsep- For Dsep 
00 rs — » 0, and for Dsep = 

'D1-D2 



<y\(J2rs 



Gi 



(Dx-D2\ „ (Dx^D2\ 



(777.2) 

rs(73i, 7)2, 0) w -/u(7)i)m(732)G23(7)i/2), 732 < 1. 
For 7)sep =0, 7)i = 7)2, rs = 1. 
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For a given AS{Di) = D\/t\ the conditional mean val- 
ues < A.Sc{D2,Dsep) > and dispersions crc{D2,Dsep) are 
described by the standard relations 



< ASc{D2,Dsep) >= 0-2rsfll/Tl, (Tc (-D2 , -Dsep ) = a2\/l - r| 

and the conditional probability of the pancake formation 
with AS{D2, Dsep) > D2/T2 under the condition AS{Di) = 
Di/ti is 



Wci> D2,T2;DuTi,D,ep)=0.5erfc 



92 
V2 



(777.3) 



92 



T2 



ri 



These relations allow us to describe the formation of pan- 
cakes with D2 < Di at the moment T2 < ti and, in partic- 
ular, to find the function < mcoid{Dsep) > and Nm{Dsep)- 
For Ti = T2 the same relations describe the accumulation of 
the smaller pancake (D2) by the larger one (Di). 

More complicated relations describe the pancake coag- 
ulation. As it was noted in Sec. 3.4 this process is charac- 
terized by the function 

(j>{qcent +D1/2) - 4>{qcent - Dl/2) 



ii>{Di,D2,Ds^p) = 



(l>{Dsep + qcent + D2/2) - (j)(D 

Sep ~\~ Qcent 

D2/2) 



D2 

and the condition of coagulation is 

IpiDl, D2, Dsep) = Dsep/r. 



{IIIA) 



To obtain the statistical description of this process the con- 
ditional mean value < ■4>'^{D-i, D2, Dsep) > and dispersion cr^ 
need to be found. Using the expressions given in Appendix 
I and notation introduced above we obtain: 
1 2 _ GoiDi) , Gq{D2) , 

l^'^lp — 7^9 1 7^9 T 



'-'2 



Go{q3i) - Go(g4i) - Go (932) + Go (942) 
D1D2 



(777.5) 



The cross correlation of functions ^ and AS{qi2), AS{q34) 
is described as following: 



(7^(72 



The dynamical characteristics of pancakes. 

A correlation of pancake velocity described by (8.4) in 
points with different transversal coordinates is convention- 
ally characterized by the correlation coefficient 

< [0(qi) - '^(q2)][0(q3) - 0(q4)] > 



y<¥(qO" 



(q2j 



>< 



(q3j 



(q4)]2 > 



Go(gi4) - Go{qia) + Go (923) - Go(g24) 
2 a/Go (ai2) Go (934) 



(7^.1) 



For small gi2,&?34 it transforms into (8.8). 

The important characteristic of the pancake is the ve- 
locity profile around the points gi&g2 under the condition 
AS{qi2) = 912/T. It is described by the mean conditional 
profile of displacement 



< Siqs) >= 



qi - 92 Gi2(gi3) - Gi2(g23) 



r 2[1 -Gi2(gi2)] 
and the conditional dispersion 

[Gi2(gi3)-Gi2(g23)]' 



?(<?3) 



1 



2[1-Gi2(gi2)] 



iIV.2) 



{IV.3) 



Using the expression (8.9) the velocity dispersion within 
pancakes can be found as follows: 



{i + zf 2 



qh 



[{l-f3{z)f+2Piz)Mq,2)- 



H^{z) " " 12 

I3'{z)f2{qi2)+T^f{z)fs{qi2)], 

^ 12 Go(?i2)-0.5gf2Gi(gi2) 



(7F.4) 



/i(gi 



9?2 



/2(gi 



/3(gi 



J12 = gi2^ 



1 - Gi2(gi2) 

3./12 
[1-Gl2(gi2)]2' 

J12 2Go(gi2) 



2[1 -Gi2(gi2)l 



qh 



dg3[Gi2(gi - g3) - Gi2(g2 - gs) 



g3iGi(g3i) - g4iGi(g4i) - g32Gi(g32) + g42Gi(g42) 



7)l7)2(T^ 

For the conditional mean value and dispersion of ip we have 

<r> _^ ( iKDi) - l^iD2)f , 2M(7)l)M(j^2) ^ 
<T^ '""y l-r? + 1 + r. )' 

<\\^_^2 Un{Di)-n{D2)f ^ 2n{D,)^{D2) \ 
cr^ J * \ 1 - rf 1+rs J ' 

Finally for the required function xiDi, D2, Dgep) we obtain 



Appendix IV 
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